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Abstract 


We study the properties of >He-*He mixture droplets using a phenomenological density functional. The structure of 
a liquid drop with and without vortices is investigated using Monte Carlo simulations. © 2000 Elsevier Science B.V. All 


rights reserved. 
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Droplets of mixtures provide a unique opportunity to 
test and enhance our understanding of finite size effects 
and the effects of boundary on the phase segregation 
between *He and “He. The study of liquid-helium mix- 
ture droplets has been a subject of great theoretical and 
experimental interest recently. In this work we present 
Monte Carlo simulations on the structure of *>He-*He 
mixture droplets with and without vortices. 

The total free energy for a liquid-helium mixture 
of *He and He can be expressed as [1] 
E = \(E3 + E, + E34)dr with E3 = (h?/2mzett)t3 (r) + 
d3|Vp3(n)? + ab3p3(r) + 2 c3ps(r)p)” + 3¢3p3(r)™™, 
E, = (h?/2m,4)|V P a(n)? + dal Vp4lr)|? + 4b3p3(r)? + 
2C4palr) p(n), E34 = d34 Vp3(r) Vpa(r) + b34p3(r)p4(r) + 
C3434 P3(r)p4(r)p(r)* where p,(r) is the particle density of 
the ‘He, p(r) is the total density p(r) = p3(r) + p4(r), and 
m; are the atomic masses. The macroscopic wave func- 
tion of *He is Y,(r)=./pa4(r)e'®. The *He density- 
dependent effective mass h?/2m3e¢¢ = h*/2m3(1 — p3(r)/ 
P3 — Palr)/p4,-)” is fitted to the experimental pressure 
dependence in the bulk liquid. In the Thomas-Fermi 
approximation, the kinetic energy densities of He can be 
expressed as 13(r) = 5 (3r?) plr)? + Te (Vps)? 
p3(r) + 3 V?°p3(r). 

From this energy functional the structure and ener- 
getics of °He-4He mixed drops without vortices have 
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been studied by means of density functional approach 
(DFA) [2]. The density distribution of *He and *He 
without vortices were obtained by solving a set of 
coupled self-consistent equations (Euler equations for 
the energy density). In all these calculations, spherical 
symmetry is assumed. In the case of dilute Bose gas 
mixture trapped in a symmetric potential, there exist 
asymmetric solutions [3]. It is not clear whether these 
symmetric solutions of *#He-*He mixed drops are really 
ground states of the system. Our calculation also pro- 
vides the density distribution at finite vorticity for the 
first time. 

Here we directly find out the lowest energy configura- 
tion with no symmetry restrictions by means of Monte 
Carlo techniques through simulated annealing. We ap- 
proximate the volume integral of the energy functional 
by a discrete sum. We sample a lattice inside a sphere of 
radius 2R consisting of 40 sites along the diameter, mak- 
ing a total of 33 398 sites. We use the Metropolis algo- 
rithm, start our simulation with high temperature, then 
decrease it to a very low temperature step by step. This is 
an annealing process and the final configuration should 
be the lowest energy configuration. 

Because of the strong interactions between the helium 
atoms, the particle number is a very sensitive function of 
the chemical potential. We thus carry out the simulations 
with fixed particle numbers. 

Fig. 1 shows the density structure of mixed droplets 
with different particle number. Because the mass of *He is 
smaller than *He, its quantum pressure is larger than that 
of *He. *He is therefore outside of *He. As can be 
seen, the structure is symmetric. The average energy 
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Fig. 1. Density distributions of *He (a,c) and *He (b, d) with 
particle numbers N3 = 1000 and N, = 3000 for (a) and (b), 
N3 = 3000 and N, = 1000 for (c) and (d). 


EN; + N4) is — 5.368 for the case of N3 = 
1000, N4 = 3000, and — 3.107 for the case of N3 = 
3000, N4 = 1000. 

The structure is different when there is a vortex in *He. 
The density distributions are given by Fig. 2. It shows 
that *He will occupy and expand the size of the vortex 
hole of liquid *He. At the same time, there are still some 
3He atoms on the outside surface of *He. This comes 
from an effective attractive potential for *He on the 
surface of *He [4,5]. With high vorticity, the vortex 
hole of He is larger and the surface attraction inside 
creates a hole in liquid *He density which is shown in 
Fig. 2(c). 
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Fig. 2. Density distributions of *He (a,c) and *He (b,d) with 
a vortex in He. The vorticity is m = 1 for (a) and (b), m = 7 for 
(c) and (d). The particle numbers are: N3 = 1000 and Ny = 3000. 
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